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Abstract 

Three types of streamline topology in a Karman vortex street flow are shown 
under the variation of spatial parameters. For the motion of dilute particles in 
the Karman vortex street flow, there exist a route of bifurcation to a chaotic 
orbit and more attractors in a bifurcation diagram for the proportion of par- 
ticle density to fluid density. Along with the increase of spatial parameters in 
the flow filed, the bifurcation process is suspended, as well as more and more 
attractors emerge. In the motion of dilute particles, a drag term and gravity 

term dominate and result in the bifurcation phenomenon. 
PACS number(s): 47.52.+j, 47.32.Cc, 05.45.-a 
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I. INTRODUCTION 



The motions of particles in a nonuniform flow have wide technological applications, such 
as to forecast chemical reactions and environmental pollution. Due to particle motion in 
the low Reynolds number category, the equation of motion for a small rigid sphere in a 
nonuniform flow field is deduced [1]. When the background flow is mainly dominated by 
large scale structures, the fluid viscosity is not included in the governing flow equation 
[2]. Related studies show that even when the background flow fields are very simple, the 
motions can have abundant phenomena. In a periodic Stuart vortex flow, depending on the 
values of parameters, the particles asymptotically concentrate along periodic, quasiperiodic 
or chaotic open trajectory [3,4]. In a cellular flow field, aerosol particles also merge into 
isolated asymptotic trajectories, which are described by slow manifolds [5,6]. Moreover, the 
method due to the Lagrange view of particles can also be applied to investigate effects of 
particle dispersion on streamwise braid vortices in a plane mixing layer [7,8]. 

In a plane wake flow behind a circular cylinder, a regular vortex street structure was 
investigated at Re = 60 — 5000. Along with the advance of experimental technique, even at 
Re = O(10 4 ), the regular vortex street is obtained by using the phase- average method [9]. 
In particular, the phenomena relating to organized vortex structure, such as, reconnection 
of vortex street [10] and emergence of three dimensional vortex structure [11], arouse wide 
interest for the transition of a plane wake flow. Recently, particle focusing in narrow bands 
near the peripheries of the vortex structures for the particle dispersion in a plane wake flow 
is observed experimentally [12]. By considering Stokes drag, the phenomenon of particle 
focusing is studied on two dimensional centre manifolds [13]. The regular Karman vortex 
street flow as a model to approach the plane wake flow and investigate the above phenomena 
plays an important role. 

For the motion of particles in regular vortex street flow, Karman vortex spacing influences 
on topological structure of background flow field. At the same time, a density ratio as a 
basic parameter may have a wide range. In this paper, we will consider streamline topology 
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and dilute particle dynamics in the Karman vortex street flow in a range of density ratio. 
In Sect. II, it is shown that for the Karman vortex street flow, there exist three types of 
global topological structure depending on the spatial parameters in flow field. Dilute particle 
dynamics in the Karman vortex street flow related to the density ratio is investigated in Sect. 
III. Effects of spatial parameters in flow field on dilute particle dispersion is determined in 
Sect. IV. Finally, a brief summary is given in Sect. V. 



II. STREAMLINE TOPOLOGY 

The stream function of Karman vortex street flow [14] is 

T / \ r , ch^(y - h/2) - cos^-x Ty , nh 

^(x, y) = —In—iJf J—L i + -%-th—, 1 

v ,yj 4vr ch^(y + h/2) + cos^x 21 I 1 K J 

where T is the strength of vortices, / and h are the streamwise and transverse spacing of 

vortices, respectively. The dimensionless quantities denoted by asterisks are introduced as 

x* = x/l, y* = y/l, h* = h/l, u* = u/U^, T* = F/iU^l) and ^* = ^/(U^l). The stream 

function (1) can be represented as 

v y ; 4tt ch27r{y* + h*/2) + cos2nx* 2 v ; 

From now on, the asterisks "*" for the dimensionless quantities in this section are omitted 

for convenience. The stream function (2) has symmetries: ^(x + 1/2, —y) = —^/(x,y), ty(x + 

l,y) = *(x,y). 

The associated velocity filed is given by 

_ <9# _ T sh2n(y-h/2) sh2n(y + h/2) I ' 

Ux ~ ~dy ~ 2^ch2n(y - h/2) - cos2nx ~ ch2n(y + h/2) + cos2nx* + 2 * ' 

Uv = — ~ — = — : r-r^ 1 — ; ; r~r^ \stn27rx, (3) 

y dx 2 l ch2Tr(y - h/2) - cos2ixx ch2n(y + h/2) + cos2vrx J v 1 

which has singularity at the vortex centers. In the case of removing the singularity by the 

Rankine vortex [15], free stagnation points in velocity field (3) consist of centers and saddle 

points, which can be respectively described by (0, |), (|, — |) and (0,—^lnc, (|, ^Inc), 



where c = 2/ shirk + shnh + ^J[2/shnh + shnh] 2 + 1. 

3 



In the symmetry ty(x + 1/2, —y) = —^>{x,y), it can be shown that there exists a zero 
streamline between two neighboring centers (0, |) and (|, — |). Under the variation of spatial 
parameter h, we take zero streamlines to investigate the evolution of global topological 
structure. The points y , corresponding to zero streamlines passing through the y coordinate, 
satisfy 

v ,y0J 4vr ch2n{y + h/2) + 1 2 v ; 

Using the Newton-Raphson and bisection methods [16], we solve Eq.(4) and plot the relation 
y ~ h in Fig. 1. 

When h — 0, the Karman vortex street returns to a periodic row of vortex pair (u x = 
T c^s2nxsh2^_ and = _r "f***^ ). In this case, a zero streamline exists in y = oo, 
as well as free stagnation points are only centers. All streamlines including an embedded zero 
streamline for h = 0.001, which approach to those for h — 0, are drawn in Fig. 2(a). Two 
neighboring vortices are divided by the zero streamline. In Fig. 1, there exists a minimum 
point (0.164,0.443434) marked by o. To understand this point, we take the derivative of 
\I/(0, y) with respect to h, as follows, 

0*(O.y) r 8 h2*(y-h/2) shMy + h/2) 9mi/rh 2 wh] n ^ 

-^h~ = -4 [ ch2n(y-h/2)-l + ch2rr(y + h/2) + 1 ~ ^ = °" (5) 

Using the Newton-Raphson and bisection methods [16], we solve Eq.(5) and plot the so- 
lutions in Fig. 3(a). Besides a trivial solution y = h/2 of Eq.(5), other relations ±y ~ h 
are symmetrical and corresponding to the maximum and minimum of ^(0, y) for h. The 
point (0.164,0.443434) exists in the curve y ~ h marked by p, so that *(0, 0.443434) is a 
maximum of ^(0, y) for h. We also plot the relation of ^(0, 0.443434) to h in Fig. 3(b). The 
maximum point q in the curve confirms to the above result. Thus, the minimum point in 
Fig. 1 corresponds to the maximum of ty(0,y) for h. 

When h increases from to 0.164, y monotonically decreases from oo to 0.443434 in 
Fig. 1. All streamlines including an embedded zero streamline for h — 0.1 are drawn in 
Fig. 2(b). In this case, saddle points emerge below or above the centers in opposition to the 



zero streamlines. The orbits passing through saddle points go around the centers and come 
back, as well as go to other saddle points, i.e., there exist homoclinic and heteroclinic orbits. 
For h = 0.164, all streamlines including an embedded zero streamline are drawn in Fig. 2(c). 
When h increases from 0.164, y monotonically increases from 0.443434 in Fig. 1. All 
streamlines including an embedded zero streamline for h = 0.3 are drawn in Fig. 2(d). When 
h increases to 0.410998, another solution of yo emerges in Fig. 1. All associated streamlines 
including embedded zero streamlines are drawn in Fig. 2(e). Zero streamlines, surrounding 
each vortex, pass through the saddle points. So, the zero streamlines are homoclinic and 
heteroclinic orbits. When h increases beyond 0.410998, the solution becomes two ones in 
Fig. 1. At the time, there exist three zero streamlines. The orbits passing through saddle 
points go around the centers and to other saddle points, but do not come back, i.e., there 
exist only heteroclinic orbits. All streamlines including three embedded zero streamlines for 
h = 0.6 are drawn in Fig. 2(f). Upper and lower vortices are divided by one zero streamline, 
as well as surrounded by others. 

In the experiment of a plane wake flow [17], it has been shown that h increases down- 
stream from 0.2 to 0.45 at Re = 140. We can thus conclude that three types of streamline 
topology can emerge in a plane wake flow. 



III. DILUTE PARTICLE DYNAMICS 



A. Governing equations 



The motion of a small spherical particle in a nonuniform flow field u is governed by the 
momentum equation [1] 

U\p P + 0.5p F )^ = f d\p P - p F )g + \<fp F ^ t + 3ndup F (u - V)f d 

+ Imv^Vf Jo* " f ) dT + f^ 3 Mu - V) x w , (6) 

where V is the velocity of the particle, d is the particle diameter, p is the density, g is the 
gravitational acceleration, v is the fluid kinematic viscosity, uj is the vorticity of the flow 



fluid, and the subscripts F and P refer to the fluid and particle, respectively. The parameter 
fd relating to Reynolds number (Rep = iHzXM) j s described [ 18,4 ^ as 

f d = 1 + 0A315Re p 82 ~° mlog ™ P , < Re P < 200. (7) 

Introducing the dimensionless quantities 5 = Pp/pf, e = 1/(0.5+5), t* — t/T,V* = V/(l/T) 
and g* = g/g (T is the particle viscous relaxation time d 2 /(18eu)), we nondimensionlize Eq. 
(6) and ignore the Basset history term. Thus, Eq. (6) can be described by 

^- = Bg + \eA 2 n . Vu + (An - V)f d + \eA(Au -V)xu. (8) 
at 2 2 

where A = U^JT/l, B — (1 — 1.5e)T 2 | and the asterisks "*" for the dimensionless quantities 

in this section are omitted for convenience. Moreover, the Reynolds number is written as 

Re P = W P /A\Au - V| (Wp = Uoodjv). 

The flow field u is chosen to be the Karman vortex street flow (3), where a parameter k 

is introduced to remove the singularities. The modified stream function (2) is described as 

T , s r , ch2ii(y — h/2) — kcos2nx Ty sh2nh . . 

\p(x v) = In — — I - (9) 

v ' yy 4tt ch2n(y + h/2) + kcos2nx 2 ch2nh + k' w 

The corresponding velocity field is given by 

_a^_r^ sh2n(y-h/2) sh2n(y + h/2) sh2nh 



dy 2 ch2ir(y — h/2) — kcos2nx ch2n(y + h/2) + kcos2irx 2 ch2iih + fc' 

Uy ~ ~~dx ~ 2^ch2ix(y - h/2) - kcos2nx + ch2n(y + h/2) + kcos27rx^ Sm2nX ' ^ 
Only when k — 1, the velocity field satisfies the Euler equation. Since the error increases 
with deviation of k from 1, we take k = 0.99 as an approximation in this simulation. 

The particle motion is described by a four-dimensional, nonlinear autonomous dynamical 
system of the form 

x = V x , 
V = V v, 

3 1 
V x = -eA 2 u • Vu x + (Au x - V x )f d + -tAuo(Au y - V y ), 

3 1 
V v = -eA 2 u • Vu y + (Au y - V y )f d - -tAuo(Au x - V x ) - B. (11) 



The parameters in Eqs. (11) can be taken as = Am/s and d = 5 x 10~ 5 m from 
[12]. Since air is chosen as the fluid media in the flow, the properties of fluid in Eqs. (11) 
are described as p = 1.225/cg/m 3 and v = 1.45 x lCT 5 m 2 /s [19]. Moreover, the parameters 
concerning the Karman vortex street flow are chosen as I = 0.1m, h — 0.3 and T = 1. 
In the following, we investigate the motion of particles under the variation of density ratio 
5. Initial values of particles are taken as points in the flow field and their corresponding 
flow filed velocity. Using a fourth-order Runge-Kutta algorithm, we integrate Eqs. (11) 
with a time size r = 0.001. After discarding transients, we plot points at x = as a 
bifurcation diagram. When the time size is changed to r = 0.01, the bifurcation diagram is 
still preserved. 

In order to understand the dilute particle dynamics, we analyze orders of magnitude 
of parameters in Eq. (8). For giving physical values in the calculation, the parameter f d 
appear to be on the order 1. When 8 is taken as the order 10 3 — 10 4 and d is fixed as 
the order 10~ 5 m, the parameters e, T, A and B appear to be on the order 10~ 3 — 10" 4 , 
10" 3 s — 10~ 2 s, 10~ 2 — 10 _1 and 10~ 4 — 10 -2 , respectively. In this case, the stress tensor term 
of fluid |eA 2 u • Vu and lift force term ^eA(Au -V)xw have smaller orders than the drag 
term (Au — V)/ d and gravity term Bg in the Eq. (8). Therefore, the Eq. (8) is dominated 
by the drag term and gravity term . 

Moreover, concerning the stability of particle orbits, we take only the drag and gravity 
terms in Eq. (8) to obtain an approximate fundamental matrix [/*. By using the straight- 
forward technique [20,21], a complete Lyapunov spectrum is determined as follows 

1 II m e° A Ut e° A Ul e2 A U* e9 II , 1 „ , , || ^.ei A E£.4 A U£.eL A [/£ 
A = lim -log - x ° \ n x ° I n x ° 3 n „ x ° 4 - = lim —Yttlog- Xl 1 Xl 2 Xl 3 



t^oo t || e° A e!j A e° A e° || n ^ nr l ~ v \\ e\ A A A 

(12) 

where A and || o || are an exterior product and a norm with respect to some Riemannian met- 
ric, respectively. After each time integration, the set of bases {ei, e^, e^, e^} is exchanged by 
using the Gram-Schmidt reorthonormalization procedure. In the calculation of the complete 
Lyapunov spectrum, the time size is taken as r = 0.001. 



B. Numerical results 



In Fig. 4, along with the increase of density ratio, a bifurcation diagram of y versus 5 is 
drawn. When 5 < 6.5 x 10 2 , the velocity of particle dispersion is very slow, so most of the 
particle trajectories are preserved near the street. When 5 > 6.5 x 10 2 , besides the particle 
dispersion, global trajectories of particles converge to two attractors: one above the street, 
the other in the street. For the attractor above the street in 5 = 6.5 x 10 2 — 4.4 x 10 3 , it 
evolves as a period-1 orbit. At 5 = 4.4 x 10 3 , the period-1 orbit bifurcates to a period-2 
orbit. In 5 = 4.4 x 10 3 - 7.0 x 10 3 , the attractor evolves as a period-2 orbit. At S — 7.0 x 10 3 , 
the period-2 orbit bifurcates to a period-4 orbit. In 5 = 7.0 x 10 3 — 7.9 x 10 3 , the attractor 
evolves as a period-4 orbit. In 5 = 7.9 x 10 3 — 1.0 x 10 4 , the period-4 orbit bifurcates to 
a period-8 orbit and further to a quasi-periodic or chaotic orbit. At 5 = 1.0 x 10 4 , a crisis 
happens, so that the quasi-periodic or chaotic orbit disappears. Moreover, at 5 = 7.3 x 10 3 , 
a period-3 orbit emerges as another attractor above the street. In 5 = 7.3 x 10 3 — 7.9 x 10 3 , 
the period-3 orbit bifurcates further to a quasi-periodic or chaotic orbit. At 5 = 7.9 x 10 3 , 
a crisis happens, so that the quasi-periodic or chaotic orbit disappears. The bifurcation 
procedure of period-3 orbit differs from that for the three dimensional Lorenz equation [22]. 
For the attractor in the street, a period-1 orbit emerges at 5 = 3.3 x 10 3 . It preserves in 
5 = 3.3 x 10 3 - 1.2 x 10 4 . 

In Fig. 5, we plot some typical examples for different values of 5. A period-1 orbit above 
the street for 5 = 1.0 x 10 3 , as well as all corresponding streamlines are drawn in Fig. 5(a). 
On the period-1 orbit, the particles move from left to right. In Fig. 6, for the period-1 orbit, 
the maximal Lyapunov exponent is -0.230, so the orbit is stable. In order to display the 
dispersion of particles in the flow field, we also plot the basin of attraction in Fig. 7(a). In 
the motion of particles, the points corresponding to those in the basin of attraction suspend 
on the period-1 orbit. At the same time, the points corresponding to those outside the 
basin of attraction escape from the central region of flow. For 5 = 3.5 x 10 3 , in Fig. 5(b), 
two period-1 orbits distribute above and in the street, respectively. On the period-1 orbit 
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above the street, the particles move from left to right. In Fig. 6, for the period-1 orbit, the 
maximal Lyapunov exponent is -0.254, so the orbit is stable. But, on the orbit in the street, 
the particles move in an opposite direction, i.e., from right to left. The maximal Lyapunov 
exponent is -1.144, so the orbit is stable. For the two period-1 orbits, the corresponding 
basins of attraction are plotted in Fig. 7(b). One is similar to that in Fig. 7(a), the other is 
distributed in two local zones. The basin for the orbit in street is surrounded by that for the 
orbit above street. For 5 = 6.0 x 10 3 , a period-2 orbit above the street and a period-1 orbit 
in the street are drawn in Fig. 5(c). On the period-2 orbit above the street, the particles 
move from left to right. In Fig. 6, for the period-2 orbit, the maximal Lyapunov exponent 
is -0.604, so the orbit is stable. But, on the orbit in the street, the particles move in an 
opposite direction, i.e., from right to left. The maximal Lyapunov exponent is -1.268, so 
the orbit is stable. The corresponding basin of attraction is plotted in Fig. 7(c), which is 
similar to Fig. 7(b). For 5 = 7.5 x 10 3 , a combined period-4 orbit with a period-3 above the 
street and a period-1 orbit in the street are drawn in Fig. 5(d). On the period-4 and period- 
3 orbits above the street, the particles move from left to right. But, on the orbit in the 
street, the particles move in an opposite direction, i.e., from right to left. In Fig. 6, for the 
period-4, period-3 and period-1 orbits, the maximal Lyapunov exponent are -1.003, -0.245 
and -1.356, respectively, so the orbits are stable. The corresponding basin of attraction is 
plotted in Fig. 7(d). The geometry of basin is different from above ones in Figs. 7(a)- (c). 
Some points near y — 1 escape from the basin of attraction in the dispersion of particles. 
For 5 = 9.5 x 10 3 , a quasi-periodic or chaotic orbit above the street and a period-1 orbit in 
the street are drawn in Fig. 5(e). On the orbit above the street, the particles move from 
left to right. In Fig. 6, for the orbit, the maximal Lyapunov exponent is 0.850, so the orbit 
is chaotic. But, on the orbit in the street, the particles move in an opposite direction, i.e., 
from right to left. The maximal Lyapunov exponent is -1.356, so the orbit is stable. The 
corresponding basin of attraction is plotted in Fig. 7(e). For the period-1 orbit in the street, 
the geometry of basin is similar to that in Fig. 7(d). But, for the chaotic orbit above the 
street, escaped points in Fig. 7(e) permeate into the basin of attraction in Fig. 7(d), so 
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that the geometry of basin is fractal. From those examples, we can conclude that along with 
the increases of density ratio, initial points distributed in the central region of flow escape 
more and more. At the same time, the particle trajectories bifurcate from periodic orbits to 
chaotic orbits. 

In order to explain the existence of attractors in the flow field, we take the period-1 and 
period-2 orbits for 5 = 6.0 x 10 3 as an example. In Fig. 8, we present the values of V x , V y , 
u y and (Au y — V y )fd — B along two orbits. For the period-2 orbit above the street, we plot 
it in two times of the streamwise periodic length in Fig. 8(a). Since V x > in x e [0, 2), the 
motional direction of particles is from left to right. (1) When a particle moves from x = 
to x = 0.486, the negative term (Au y — V y )f d — B causes decrease of V y from 0.050 to -0.285. 
At the same time, it leads to increase the term (Au y — V y )f d — B. (2) When the particle 
moves from x = 0.486 to x — 0.930, the positive term (Au y — V y )f d — B causes increase of V y 
from -0.285 to 0.404. At the same time, it leads to decrease the term (Au y — V y )fd — B. (3) 
When the particle moves from x = 0.930 to x — 1.554, the negative term (Au y — V y )f d — B 
causes decrease of V y from 0.404 to -0.115. At the same time, it leads to increase the term 
(Au y — V y )f d — B. (4) When the particle moves from x = 1.554 to x — 1.946, the positive 
term (Au y — V y )f d — B causes increase of V y from -0.115 to 0.067. At the same time, it leads 
to decrease the term (Au y — V y )f d — B. (5) When the particle moves from x = 1.946 to 
x = 2, the negative term (Au y — V y )f d — B causes decrease of V y from 0.067 to 0.050. At 
the same time, it leads to increase the term (Au y — V y )f d — B. According to the periodic 
boundary condition, when the particle reaches x — 2, it goes back to x — 0. Therefore, 
the combination of the drag term and gravity term in the vertical direction has a periodic 
vibration along with the variation of the vertical velocity of particles and makes the period-2 
orbit. In Fig. 8(b), for the period-1 orbit in the street, we plot it in the streamwise periodic 
length. Since V x < in x G [0, 1), the motional direction of particles is from right to left. 
In a similar way to Fig. 8(a), the term (Au y — V y )f d — B brings into a periodic vibration 
of particles and makes the period-1 orbit. In Ref. [13], by considering the drag term, the 
essential dynamics takes place on the two-dimensional centre manifolds. To compare with 
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the result, we find the bifurcation process disappears when eliminate the gravity term in 
Eq. (8). From the above observation, we can conclude that the drag term and gravity term 
lead to the bifurcation behavior in dilute particle dispersion. The motional direction of 
particles and distribution of u y determine the vertical position of attractors. 

IV. EFFECTS OF SPATIAL PARAMETERS IN FLOW FIELD ON DILUTE 

PARTICLE DISPERSION 

In the Sect. II, we show three types of streamline topology in the Karman vortex street 
flow at different spatial parameter h. In the Sect. Ill, we present the dilute particle dynamics 
for the flow field with h = 0.3. In the following, we increase the spatial parameter h to 
consider its effects on dilute particle dispersion. 

For h = 0.410998, along with the increase of density ratio, a bifurcation diagram y 
versus 5 is drawn in Fig. 9. The global evolution is similar to that in Fig. 4. However, the 
bifurcation of period-1 attractor above the street is delayed. The period-1 orbit above the 
street bifurcates to a period-2 orbit at 5 = 6.6 x 10 3 . The period-2 orbit bifurcates to a 
period-4 orbit at 5 = 1.04 x 10 4 . In 5 = 1.04 x 10 4 — 1.2 x 10 4 , the period-4 orbit bifurcates 
further to a quasi-periodic or chaotic orbit. A period-3 orbit above the street emerges at 
5 = 7.0 x 10 3 and bifurcates in 5 = 7.0 x 10 3 - 8.3 x 10 3 . At 5 = 8.3 x 10 3 , a crises 
happens, so that the quasi-periodic or chaotic orbit above the street disappears. Moreover, 
the occurrence of period-1 attractor in the street is shifted earlier at 5 = 2.0 x 10 3 . A 
bifurcation of period-3 orbit in the street emerges in 5 = 9.4 x 10 3 — 1.1 x 10 4 . 

For h = 0.6, along with the increase of density ratio, a bifurcation diagram y versus 
5 is drawn in Fig. 10. The global evolution is similar to that in Fig. 9. However, the 
bifurcation of period-1 attractor above the street is delayed. The period-1 orbit above 
the street bifurcate to a period-2 orbit at 5 = 9.1 x 10 3 . The period-2 orbit is preserved in 
5 = 9.1 x 10 3 — 1.2 x 10 4 . At the same time, the occurrence of period-1 attractor in the street 
is shifted earlier at 5 = 2.0 x 10 2 . Moreover, at 5 = 7.6 x 10 3 , a period-4 and a period-3 orbits 
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emerge above the street. They bifurcate further to two quasi-periodic or chaotic orbits and 
disappear at 5 — 8.4x 10 3 and 5 = 8.5x 10 3 , respectively. In 5 — 8. Ox 10 3 — 9.4x 10 3 , a similar 
bifurcation process of more attractors emerges in the street. In 5 — 8.0 x 10 3 — 8.8 x 10 3 , a 
period-4 orbit preserves in the street. In S = 8.4 x 10 3 — 9.4 x 10 3 , a bifurcation of period-3 
orbit exists in the street. 

Thus, along with the increase of h, a period- 1 orbit above the street bifurcates slower 
to a quasi-periodic or chaotic orbit, as well as a period-1 orbit in the street emerges earlier. 
Besides the main attractors in the flow field, more and more local attractors also appear 
above and in the street. 

V. CONCLUSION 

Under the variation of spatial parameters, we have shown three types of streamline 
topology in a Karman vortex street flow. For the motion of dilute particles in the Karman 
vortex street flow, there exist a route of bifurcation to a chaotic orbit and more attractors 
in a bifurcation diagram for the proportion of particle density to fluid density. Along with 
the increase of spatial parameters in the flow filed, the bifurcation process is suspended, as 
well as more and more attractors emerge. In the motion of dilute particles, a drag term and 
gravity term dominate and result in the bifurcation phenomenon. 
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FIGURES 



Fig. 1 Points of zero streamlines in the y coordinate. 

Fig. 2 All streamlines including corresponding embedded zero streamlines with centers 
marked by dots and saddle points marked by crosses for (a) h = 0.001; (b) h = 0.1;(c) 
h = 0.164; (d) h = 0.3; (e) h = 0.410998; (f) h = 0.6. 

Fig. 3 (a)A relation y ~ h for f mai / m i n (0, y); (b) Stream function relating to h at 
(0,0.443434). 

Fig. 4 A bifurcation diagram with h = 0.3 for a continuous range of 5 showing the 
vertical position of particles at x — 0. The dense bifurcation zone 5 G [7000, 10000] above 
the street is enlarged and redrawn in the bottom-right corner of the figure. 

Fig. 5 Typical trajectories in Fig. 4 for (a) 5 = 1.0 x 10 3 and all streamlines; (b) 
5 = 3.5 x 10 3 ; (c) 5 = 6.0 x 10 3 ; (d) 5 = 7.5 x 10 3 ; (e) 5 = 9.5 x 10 3 . 

Fig. 6 Variation of the maximal Lyapunov exponent \ ma x with 5. Besides the period-3 
orbit for 5 = 7.5 x 10 3 labeled by a triangle, the attractors above and in the street are 
marked by dots and crosses, respectively. 

Fig. 7 Basins of attraction corresponding to the typical trajectories in Fig. 5 for (a) 
5 = 1.0 x 10 3 ; (b) 5 = 3.5 x 10 3 ; (c) 5 = 6.0 x 10 3 ; (d) 5 = 7.5 x 10 3 ; (e) 5 = 9.5 x 10 3 . 
In Figs. 7(b)-(e), besides two local zones marked by thinner lines for the basin of period-1 
orbit in the street, the global zone presents the basin for attractors above the street. 

Fig. 8 Distribution of V x , V y , u y and (Au y — V y )fd — B along x on (a) a period-2 orbit 
above the street; (b) a period-1 orbit in the street for 5 = 6.0 x 10 3 . 
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Fig. 9 A bifurcation diagram with h = 0.410998 for a continuous range of 5 showing the 
vertical position of particles at x — 0. The dense bifurcation zone 5 G [7000, 8300] above the 
street is enlarged and redrawn in the bottom-right corner of the figure. 

Fig. 10 A bifurcation diagram with h = 0.6 for a continuous range of 5 showing the 
vertical position of particles at x — 0. The dense bifurcation zone 5 G [7600, 8500] above the 
street is enlarged and redrawn in the bottom-right corner of the figure. 
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